Highly Dispersive Scattering From Defects In Non- Collinear Magnets 
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We demonstrate that point-like defects in non-collinear magnets give rise to a highly dispersive 
structure in the magnon scattering, violating a standard paradigm of its momentum independence. 
For a single impurity spin coupled to a prototypical non-collinear antiferromagnet, we find that the 
resolvent is dominated by a distinct dispersive structure with its momentum-dependence set by the 
magnon dispersion and shifted by the ordering vector. This feature is a consequence of umklapp 
scattering off the impurity-induced spin texture, which arises due to the non-collinear ground state 
of the host system. Detailed results for the staggered and uniform magnetization of this texture as 
well as the T-matrix from numerical linear spin- wave theory are presented. 
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Introduction. — Electron localization [lj, paramagnetic 
impurities in superconductors [2j, and the orthogonal- 
ity catastrophe [3j, all attest to the fundamental im- 
portance of impurities as probes of quantum many-body 
systems. Major research effort in cuprate superconduc- 
tors has led to extensive studies of impurities in the 
square-lattice Heisenberg antiferromagnets (HAFs), un- 
covering new universality classes for disorder-driven tran- 
sitions @HS], impurity-induced magnetic order [9], frac- 
tional Curie response [TQl [TTj, and anomalous low-energy 
magnon scattering [12} [T3j . 

While the square-lattice HAF is unfrustrated and has 
a collinear ground state, defects in non-collinear and frus- 
trated quantum magnets have come into focus only re- 
cently, displaying an even richer physics. This includes 
frustration release, dimer freezing, and mutual impu- 
rity repulsion [T4HT6] , valence bond glass states [T71 [T8] , 
emergent gauge-flux pinning [19J, breakdown of linear re- 
sponse [20J, fractional impurity moments, and — the pri- 
mary topic of this Letter — spin textures [2TH23] . 

Impurity-induced spin textures are a genuine hallmark 
of non-collinear magnetic order and can be understood 
on a purely classical level. Removing a spin from the 
host, or adding an extra defect spin, locally perturbs the 
balance of exchange fields and requires the surrounding 
spins of the non-collinear host to readjust their directions 
recursively, resulting in a long-ranged modification of the 
canting angles, i.e., a texture [2TH23] . A ID sketch of this 
is shown in Figs, [ljb) and (c) for the field-induced non- 
collinear state coupled to an impurity spin. The readjust- 
ment effect is absent for collinear order, where impurity 
spin simply co- aligns with the host, as in Fig. [TJa). In 
contrast to that, the texture implies a fractional screen- 
ing of the impurity moment |22j . The real-space decay 
of the texture depends on the nature of the non-collinear 
state. In a field-induced canted states, textures decay ex- 
ponentially on a length scale inversely proportional to the 
external field [21J. In frustration-induced non-collinear 
states, Goldstone modes lead to an algebraic decay of 
the texture [22H24) . 



In this Letter we advance the field beyond previous 
studies, which have focused on the static properties of 
defects, and investigate magnon impurity-scattering in 
non-collinear magnets. To be specific, we consider the 
field-induced canted state of the square-lattice HAF with 
an additional defect, namely an extra out-of-plane spin 
interacting by an exchange coupling with one of the host 
spins. We discover a phenomenon rather surprising, if 
confronted with conventional expectations for the scat- 
tering amplitude from a point defect, which is either 
momentum-independent altogether, aside from the triv- 
ial transformation of the excitation basis, or contains only 
a broad momentum modulation due superposition of a 
few partial waves. Instead, the scattering amplitude dis- 
plays a strongly dispersive feature, clearly tracing the 
magnon dispersion shifted by the magnetic ordering vec- 
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FIG. 1. (color online) (a) Impurity spin coupled to a collinear 
state: all spins co-aligned, (b) Homogeneous canted state 
in external field h z . (c) Impurity spin coupled to the canted 
state: host spins readjust, creating a texture, (d) ID sketch of 
umklapp scattering by the texture, which generates staggered 
^-component of the effective field with the wave vector Q = 
7r. (e) Solid black line: magnon dispersion; blurred red line: 
dispersive peak in scattering amplitude. 
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tor. We show that this effect is an unequivocal conse- 
quence of the spin texture. Intuitively, an effective stag- 
gering of the magnetic field is generated by the texture, 
made explicit in Fig. [TJd). This serves as a potential for 
umklapp scattering of magnons, which, in turn, leads to 
the central new feature in the T-matrix — a momentum- 
dependent resonance. In the following, we provide the 
detailed arguments for this result, which should remain 
valid for a wide class of frustrated non-collinear systems, 
and suggest experiments to test this prediction. 

Model — We consider the square-lattice HAF at T = 
in an external field, coupled to an impurity spin S' 
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where (Im) are the nearest-neighbor bonds of the square 
lattice, the exchange couplings of the host (Jo) and host- 
to-impurity (J) are antiferromagnetic. The gyromagnetic 
ratio is identical for all spins and is included into the 
magnetic field h. In the following, we set Jo = 1. 

The spin configuration that minimizes the classical en- 
ergy of model at h^O corresponds to an inhomoge- 
neous distribution of spin tilt angles Oi out of the xy- 
plane where ordering occurs at h = 0, see Fig. [1] For a 
1/S expansion, we align the local spin quantization axis 
on each site in the direction given by the local canted 
frame [25} [26j. The rotation of spin components from 
the laboratory frame (xo, yo,zo) is given by Sf° =Sf and 
s xo(z ) =s x(z) sin0l ±S^ x) e^ cos ft, where Q = (tt, tt) 
is the Neel ordering wave- vector. The transformation is 
the same for the impurity spin as it can be seen as a 
neighbor of the site / = 0, which it is coupled to. 

Expressing the spin operators in terms of Holstein- 
Primakoff bosons, Hamiltonian is transformed into 
a series H = H c \ass + T~Li + % + ••• with decreasing pow- 
ers of S(S') and increasing number of boson operators. 
Each term in this series depends on all 6^ and H c iass is 
the classical energy [27J. The harmonic spin- wave term is 
H.2 and stability requires 'H 1 =0. Equivalently, the ground 
state must minimize 'Hciass, i-e. 3% c i ass /<90/n = 0. With- 
out the impurity, all Oi = sin -1 (h/h s ) with the saturation 
field h s = SS [25J. With the impurity, minimization gives 
a set of nonlinear coupled equations, which determine the 
inhomogeneous distribution of the local tilt angles Oi — 
referred to as the texture hereafter. 

In what follows, we study the properties of this tex- 
ture numerically in finite N x N clusters with periodic 
boundary conditions. First, we briefly address its static 
properties and then turn to its quantum dynamics using 
numerical real-space diagonalization of H2- 

Classical texture. — The spatial extent and field- 
dependence of the texture can be described in terms of 
the staggered z-component of the magnetization rn z tagri 
obtained from the set of sin(ft). Our results, inset (a) 
of Fig. [5] and [27J, largely corroborate earlier findings 
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FIG. 2. (color online) Impurity magnetization m tmp vs h for 
J=l, 3, and 6 in N = 72 cluster and for J=l in N = 64 
cluster (dashed). Insets: (a) Local magnetization Amf = 
(sin(0o) — sin(0z)) in a 21x21 section of the N — 12 cluster, 
for J=l, h=0A. (b) Local magnetization mf at the distance 
(A/2, A/2) from / = in A = 8 and 16 clusters for J=l. 



of [21], where m z stag Yi was investigated by a continuum 
theory and quantum Monte Carlo for a different impurity 
type. In particular, the texture decays exponentially at 
|r;|^>l, consistent with the impurity not coupled to the 
Goldstone mode of the host system. 

Fig. [2] shows another characteristics of the texture: 
the impurity contribution to the uniform magnetization 
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host 



vs field for several values of the cou- 
pling J. Here m z =^2 l Sf is the uniform magnetization 
including S[ z and vn z host = $i 18 ^ na ^ °f tne nost 
in the absence of impurity. We use S=S'=1 hereafter. 
Defining the impurity susceptibility as Ximp = dmi m p/dh, 
Fig. [2] shows several regimes of screening of the impurity 
by the texture: partial, complete, and overscreening, as 
evidenced by Ximp > 0, ^0, and < 0, respectively. This 
is consistent with a field-dependent fractional effective 
impurity spin [22], and is in a stark contrast with the 
collinear HAFs, where classical mi mp = S' . The impurity 
magnetization is critical at h s because the susceptibility 
of the host diverges as ~ In \ h — h s \ [25J. Fig. [2] also shows 
that the saturation in the system with impurity occurs 
above h s of the pure host and that finite-size effects are 
negligible for the clusters and field ranges that we use. 

For completeness, we note that the impurity-induced 
classical texture behaves singularly at h —> 0, although 
in a field range of measure zero in the thermodynamic 
limit — an effect also noted in [20) [28] . In a finite 
system, the energy gain of the canted state in Fig. [T] 
AE - -N 2 h 2 /(8S), is less than that of the state in 
which the Neel order of the host and the impurity spin 
both fully align with the field, AE — —hS'. Thus, at 
h = 0+ host spins are aligned (anti-aligned) with the field, 
Sf = ±5. A spin-flop crossover to the textured state oc- 
curs at h c - 8SS f /N 2 H> as A^oc. Inset (b) of Fig. [2] 
displays this behavior on judiciously small systems by 
monitoring the magnetization mf of a spin at the largest 
geometrical distance from the impurity. 
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T-matrix. — We now turn to the spectral properties 
of the system. Because the texture breaks trans- 
lational invariance, the Bogolyubov transformation of 
H2 has to be performed numerically [29J. The para- 
unitary, 2(7V 2 + l)x2(7V 2 + l) matrix U of this transfor- 
mation maps the local Holstein-Primakoff bosons = 
[a[ , . . . , a^ N2 , a\ , a t , . . . , a N2 , a i ] onto Bogolyubov bosons 
= . . . , 5] v - 2+1 , b 1? . . . , 6jy2 +1 ]j whose Hamiltonian, 
H = b t Eb/2, is diagonal. The eigenenergies E n are 
all positive except for Ej = of the Goldstone mode. 
The Green's function in the b-basis is also a diagonal 
2(A/' 2 +l)x2(A/' 2 +l) matrix G 5 (» = [2P-E]- 1 , where P 
is the para- unit matrix with 1(— 1) in the upper (lower) 
half of its diagonal. The Green's function of the local 
Holstein-Primakoff bosons is G a (z) = (Ut)- 1 G 5 (z)U- 1 . 

However, to formulate the scattering problem for the 
impurity-induced texture, the proper basis is that of the 
Bogolyubov magnons of the uniform host, which describe 
the incident and scattered magnons as plane-wave eigen- 
states of momentum k. Thus, we first Fourier transform 
the matrix elements of G a (z) of the local host bosons 
to k-space, yielding a matrix G k , k (z). Second, the host 
boson terms of this matrix are mapped onto the basis 
of the Bogolyubov magnons b^ = [b k ,6 k ] of the uniform 
host, using the known parameters of the transformation, 
Uk and for the square-lattice HAF in a field [26) [27]. 
This yields a matrix Green's function with three 2x2 
substructures made from blocks of rank N 2 x A/" 2 , 1, and 
N. They correspond to the dressed (i) host magnon, 
(ii) impurity, (iii) and magnon-impurity Green's func- 
tions Gk'k(4 Gi(z), and G k ^(z), respectively. 

Altogether, starting from the numerical solution of the 
classical texture, followed by the real-space diagonaliza- 
tion of the harmonic problem, and Bogolyubov trans- 
formation onto the uniform host, we obtain the dressed 
magnon Green's function Gk'k(^). On the other hand, 
Gk'k(^) can be written in the conventional form 

G k / k (*) = S^Gi(z) + Gl,(z)Tm(z)Gl(z) , (2) 

where G k (z) is the diagonal 2x2 Green's function of 
the uniform host magnons with G k (2) = G k (— 2) = 
[z — £k] -1 and £k is the magnon energy. Using Eq. ([2]), 
we can now extract the scattering matrix Tk'k(^) from 
Gk'k(^)- For the remainder of this work we focus on the 
diagonal elements of the T-matrix, Tkk(^) = Tk'k(^)4'k 

T kk (^) = [Gi(z)}- 2 G kk (z) - [Gi(z)}- 1 , (3) 

which suffice to state our main findings. 

No texture test. — First, we demonstrate the feasibility 
of obtaining the T-matrix from Eq. ([3| numerically. For 
that purpose, we solve a complementary artificial prob- 
lem, in which we neglect the feedback of the impurity on 
the host spins, i.e., spins in the plane retain their homo- 
geneous field-induced canting of Fig. [TJb) and no texture 
is created. While such a reference state is, of course, 




FIG. 3. (color online) Analytical and numerical results for the 
T-matrix spectrum in no-texture case, J=l, h—1. Homoge- 
neous canting angles of the host spins 6~ 0.1253. Impurity 
canting angle 0i~ 0.8729 as in the actual texture. Thin solid 
blue: exact —Im t 11 (z=lj) [27]. Impurity resonance energy 
at e ~ 1.308, anti-bound state at ~ 4.246, van Hove singu- 
larities at ujgap — 1 and oj max & 4 are clear. Thick red solid: 
numerical — Im£ k } =k==0 (2=u;-|-i0. 05) for A=64. Inset: numer- 
ical — Im t^ k (z= ^+^0.05) along the k-path of Fig. [4] 

unstable as H c iass is not at its minimum, it permits an 
analytical solution of the scattering problem of H2, de- 
tails of which are provided in [27J. The analytical solution 
can be compared to the T-matrix obtained from the nu- 
merical procedure described above. In the following we 
consider the resolvent, i.e., the T-matrix stripped from 
the matrices of the Bogolyubov basis transformation 

t k / k (*) = (Bj^^TVk^XBk)- 1 , (4) 

where £ k 1(22) =w k and 5 k 2(21) =v k [27J. 

The analytical result for the the resolvent spectrum, 
— Imt k } k (z), is plotted in Fig. [3] vs frequency uj. Nat- 
urally, tk'kW = t(z) is momentum independent |27| . 
This is an expected behavior for scattering from point- 
like defects and is similar to scattering from vacancies in 
collinear H AFs [12j [13] , where the resolvent shows some 
broad k-modulation from superposition of a small num- 
ber of partial waves. The inset of Fig. [3| shows — Im t\^(z) 
obtained numerically from ^ and ([4]) along the path in 
k-space shown in Fig. [4j Clearly, it is also momentum 
independent. In addition, analytical and numerical re- 
sults, if evaluated on the finite clusters of the same size, 
agree to within numerical precision |27| . 

Finally, Fig. [3] demonstrates the spectral resolution we 
can obtain from the numerical procedure in an N = 64 
cluster with a minimally acceptable imaginary broaden- 
ing. One can see, that the numerical scattering ampli- 
tude has all the features of the analytical one: the im- 
purity resonance, the shallow spin- wave continuum, and 
the anti-bound state above the upper edge of the spec- 
trum [30 J. Fine details, such as the anti-bound state gap 
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FIG. 4. (color online) The T-matrix resolvent spectra 
— Imt^k(^— <^+i0.05) in N = 64 cluster vs k and cj, J = 1, 
/i=l, 2, and 4 along the depicted k-path. Panel (a): contour 
plots superimposed with the shifted magnon dispersion £k+Q 
(red-yellow dots). Panel (b): 3D plots. 



and the non-analytic van Hove singularities are smeared 
out. Improving this with systems sizes beyond N = 70 
is impractical because of the large memory requirements 
for the non-sparse 2(7V 2 + 1) x2(7V 2 + l) matrices. 

Dispersive resonance. — We now consider the scatter- 
ing T-matrix for the true ground state of the system with 
the spin texture. An analytical solution is not possible 
in this case. With the feasibility of the numerical pro- 
cedure established, we evaluate the T-matrix from Eq. 
(|3| using 0{iy from the minimization of H c iass as an in- 
put to the Bogolyubov transformation. Representative 
results are shown in Fig. [4j Removing the k-dependence 
due to transformation of the basis from Tkk(^) as in Q, 
we show — Im£ kk (z) as a function of uo and k along a 
high-symmetry path in the Brillouin zone and for several 
values of the magnetic field. 

In a sharp contrast to the no-texture case, Imt kk (z) 
reveals a clear dispersive feature. The localized impurity 
resonance in Fig. [3] is now visible only as a faint maxi- 
mum and is completely overshadowed by the dispersive 
resonance. Such a result is completely unexpected for 
the point-like impurity coupled to the Heisenberg model 
0. Direct comparison in Fig. [4^a) shows that the k- 
dependence of the dispersive resonance closely follows the 
spinwave dispersion £k+Q, folded by the ordering vector 
Q = (7r, 7r). As one can see, the resonance is most sharply 



defined for small fields and gets washed out at higher 
fields. We find the dispersive feature to be prominent 
regardless of the system size or the impurity coupling J. 

It is reasonable to suggest that the dispersive resonance 
is a natural outcome of the scattering from an extended 
region of the impurity-induced texture, arising due to 
non-collinearity of the state. This can be understood 
qualitatively from Fig. [TJc), which shows that the impu- 
rity spin has a component that acts as a local field in the 
direction perpendicular to the homogeneous field-induced 
canting. Because of that, the spins of Fig. [TJb) are per- 
turbed from their local reference frames by the staggered 
transverse effective field. Then the spin-wave part of the 
Hamiltonian can be written as % = ^H h ' 1 -\-7i stag , where 
q_jh,i con tains the homogeneous canting of spins and the 
point-like impurity scattering as in the no-texture case, 
while / H sta 9 is inhomogeneous with staggered matrix ele- 
ments, which decay on the length scale set by the texture. 

Because of the staggering, magnons must experience 
an umklapp scattering potential that can be approxi- 
mated, for an extended region of the texture, as / H sta Q ~ 
J^ k Wk^ k +Q^k- Here, a qualitative analogy can be 
drawn with the ID Kronig- Penney model whose T-matrix 
is dispersive and has a pole close to the zone-folded band 
Ek+Q [31J. Because of the finite spatial extent of the tex- 
ture, the dispersive resonance must be broadened. This 
is consistent with the increase of the broadening in Fig. [4] 
at higher fields where the size of the texture shrinks. This 
may imply a nontrivial behavior of the T-matrix in the 
limit of h — >• where the texture becomes quasi-long- 
ranged as in magnets that are non-collinear in zero field. 
We note that the impurity scattering does not lead to 
overdamping of the Goldstone mode, i.e., the spectral 
density at low energies in Fig. [4] does not occur at the 
ordering vector Q. 

Our results are of a direct relevance to the excitation 
spectra of non-collinear magnets with a low concentration 
x of naturally occurring or deliberately doped impuri- 
ties. Since the magnon self-energy is simply proportional 
to the diagonal element of the T-matrix via Sk(cj) ~ 
xTkk(^), one may expect to observe an anomalous k- 
dependent broadening of the spectrum where over- 
lap with £k+Q and an equally unusual field-dependence 
of such a broadening. These and other features should 
be observable by inelastic neutron scattering and specific 
predictions will be subject of future work. 

Conclusions. — To conclude, we have presented strong 
evidence for a highly anomalous static and dynamic re- 
sponse of non-collinear antiferromagnets to doping by 
point-like defects. The scattering amplitude exhibits fea- 
tures that are strikingly different from the usual s-wave 
scattering and include a highly dispersive resonance due 
to an impurity-induced texture. This result should be 
valid for the broad class of non-collinear magnets. Fur- 
ther theoretical and experimental studies seem highly de- 
sirable. 
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Classical Hamiltonian 

The classical energy of the square-lattice Heisenberg 
AF in a field with an out-of-plane impurity, model (1) of 
the main text, is 

^class = ~S 2 Jo ^ COS (^ +Qm) ~ sin S l 

(Im) I 

-JSS' cos(6> + - hS' sin 0* , (5) 

where (Im) denotes bonds. The impurity site i is coupled 
to site I = of the host and # n 's are the tilt angles out 
of the x?/-plane, see Fig. [5] The classical ground state 
is obtained by numerical minimization of H c iass? i-e. by 
solving the set of equations cW c i ass / dO ^ = 0. The re- 
sulting spin configuration corresponds to an inhomoge- 
neous distribution of spin canting, parametrized by the 
local tilt angles 0j, i.e. the texture. Our Fig. [6] exhibits 
one of the quantities that can be used to analyze the spa- 
tial extent and other characteristics of the texture: the 



staggered component of the magnetization m 



stag,r' 
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main panel displays m z stag r along the x-axis, at location 
r = (r, 0) off the impurity, where 
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(6) 



for two different values of impurity coupling J and for 
various field strengths. The impurity is coupled to site 
R = (0, 0) of a cluster with N = 70. Clearly, the spa- 
tial extent of the texture increases as the field h —> 0. 
Inset (a) demonstrates that the texture decays exponen- 
tially [lj at r ^> 1. This behavior is expected because 
the impurity is not coupled to the Goldstone mode of 
the host. This finding is consistent with earlier work [2 J 
where m z stag Y was investigated for a different type of im- 
purity (vacancy) by a continuum theory and quantum 
Monte Carlo. We do not observe exact exponential be- 
havior in the low- field regime, most likely because of the 
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FIG. 6. (color online) The staggered component of mag- 
netization m,stag,r ^ vs distance r from the origin for two 
different values of impurity coupling J and for various field 
strengths. The impurity is coupled to site R = (0, 0) of a 
cluster with N = 70. Inset (a) same on the semi-log plot. 
Inset (b) m z stag ^/h vs r • h. 

finite cluster size. By varying h we find that Tn z stag r /h 
scales almost perfectly with r • h, see inset (b) of Fig. [6| 
again in agreement with [2 J . One potential reason for the 
visible deviation from scaling is the field dependence of 
the transverse susceptibility x± an d the spin stiffness p s , 
neglected in the continuum description of [2]. 



Harmonic part of the Hamiltonian 

Within the 1/S expansion, the local spin quantization 
axes on each site are aligned in the direction given by 
the local canted frame with the angle 0j, obtained from 
the minimization of the classical energy. The subsequent 
Holstein-Primakoff bosonization of spin operators yields 
the harmonic spin- wave Hamiltonian %2- The term lin- 
ear in bosonic operators, Hi, vanishes identically upon 
minimization of H c iass- 

The spin-wave Hamiltonian of the host reads 
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FIG. 5. (color online) A sketch of the spin configuration. 



where the summation is over the lattice sites I and the 



2 



nearest neighbors 5 and the shorthand notation 9i$ = 
0i + Qi+s has been introduced. The impurity part of the 
Hamiltonian is 



l/imp 

n 2 



J cos 0q 7 



sin 



( a o a < + H - c -) 
2" ( a o a * 



h sin Oj 



(8) 



H.c. 



where #oi = 9q + Oi. The first line contains a potential- 
like energy shift for the magnon on the site coupled to 
the impurity (I = 0) and the local energy of the impu- 
rity magnon, while the rest of the Hamiltonian describes 
various transitions between the two. 

In the textured case, analytical diagonalization of the 
Hamiltonian ([7| and ^ is not feasible, so we perform the 
Bogolyubov transformation numerically. Subsequently, 
the T-matrix is extracted from the Green's function writ- 
ten in the basis of Bogolyubov magnons of the uniform 
system (no impurity), as described in the main text. 



No texture test 

The feasibility of extracting the T-matrix from the full 
Green's function numerically can be demonstrated for a 
complementary artificial problem, for which an analyt- 
ical solution can also be found, and by comparing the 
results of the two approaches. Here we formulate such a 
problem by neglecting the feedback of the impurity onto 
the host spins, so that no texture is created. While such 
a reference state is unstable as rl c \ ass is n °t minimal, it 
permits an analytical solution of the scattering problem. 



Hamiltonian 

In the no-texture case, the canting of the host spins 
is uniform, 9\ = 9, with the canting angle found from 
the energy minimization of the system without impurity: 
sin# = h/h s with the saturation field h s = 8SJo- The 
subsequent diagonalization of the host Hamiltonian ^ 
is straightforward [3j |4] and leads to 



^host 



where = 4SJo + 7k) (1 — 7k cos 29) is the spin- 
cos k y ). The Bo- 



wave dispersion and 7k = ^(cos/c^ 
golyubov transformation from 
and is written as 



ikr; 



B k 1< 



aj to 6 k 



(10) 



where 



Bl 



Vk «k 



(11) 



is a convenient 2x2 notation. The u-v factors are: 
^k = VT^k + £k)/2e k , ^k = sign(7 k ) ^(Ak - £ k )/2£ k , 
with A k = 4SJo(i + 7k sin 2 9) 0]. 

Using these notations and that Bk = B^, the impurity 
part of the spin- wave Hamiltonian |8j) can be written as 



ea]a- 



k,k' 



with $ ? 



T E *lc'Bk'B k * k + ]T ^BkV^ , 

(12) 



Vi = V 1 



f 9 
9 f 



where potential- like scattering amplitude Vq = JS cos #oi? 
magnon-impurity transfer amplitude V\ = Jy/SS' \ ma- 
trix elements / = sin 2 (#oi/2) and g = — cos 2 (#oi/2), and 
the angle 9oi = 9 + 9i as before. With the notations (11), 



the last term of the magnon-impurity scattering in (12) 



implicitly contains its own conjugate. The impurity en- 
ergy is £ = JS r cos 9oi + h sin 9{ with an impurity canting 
angle 9{ which is a free parameter. In the following, we 
fix the latter to its value obtained numerically for the 
problem with the texture. One might also chose 9i to 
minimize the energy of the impurity spin in ([5| while 
keeping the host canting angle homogeneous. This yields 
tan#i = (8 Jo/ J— 1) tan 0, which is numerically close to 
the problem with the texture. 



Green's function and T-matrix 

The imaginary time Green's function of the 
magnons is a 2x2 matrix which can be written 
as a direct product G k / k (r) = -(T T (* k/ (r) <g> 
using ( pT| ). Its Fourier transform is G k 'k(^n) = 
exp(icj n r)Gk'k(T)^T, where iuo n = i2nnT, which we 
replace by a complex variable z hereafter. The Green's 
function in the presence of impurity scattering can be 
expressed through the T-matrix 

G Vk (z) = <5 k , k G k (^) + G k ,(z)T k , k (z)G k (^) , (13) 

where the noninteracting Green's function G^(z) is 

1 



Gi(z) = 







(14) 



(9) Using the impurity scattering terms in ( 12 ), the T-matrix 



follows from an infinite sequence of the two-component 
vertex function, one component from the potential-like 
scattering Fo-term and the other from the magnon- 
impurity scattering Vi-term. Figs.[7^a) and (b) show the 
vertex and the sequence, respectively. From the structure 
of the vertex 



r k / k (^) = Bj c/ r(^)B k , 

where the k-independent vertex is 

r(z) = v i + v 1 G° i (z)v 1 , 



(15) 



(16) 
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FIG. 7. ( colo r online) (a) Composite vertex of impurity scat- 
tering in ( |12| ). (b) T-matrix sequence with that vertex. Dot- 
ted and solid lines are the noninteracting impurity resonance 
and the host magnon Green's functions, respectively. 



and G®(z) is the local Green's function of the impurity 

state, 



G?(*) 





i 

' Z + £ 



(17) 



it is straightforward to see that the T-matrix has no In- 
dependence aside from the trivial basis change Bk and 
Bk' to Bogolyubov bosons 



T k / k (z) =Bkt(z)B k . 



(18) 



This is a standard feature, not only for the particular 
type of impurity in antiferromagnets that we investigate, 
but for the point-like scatterers in general. Stated differ- 
ently, the actual scattering resolvent t(z) is completely 
momentum independent Using (15), ( [l6| ), and Fig. (7^b) 
one readily arrives at the final expression for the resolvent 



t(z) = 

where the matrix J?(z) 
F(z) = 



T(z) 



■F(z)T(z)> 



G (z) F (z) 
F (z) Go(-z) 



(19) 



is built from the local Holstein-Primakoff Green's func- 
tions of the host 



Go(z) = J2 

k 

Fo(z) = £ 



Vk^k 



z + £k 
iik^k 



£k 



£k 



(21) 



This concludes the formal solution of the scattering prob- 
lem. The remaining step is the momentum integration in 
Gq(z) and Fq(z). In the next section, we summarize our 
analytical results for this integration. 



Local Green's functions 

Since in the scattering problem we are interested in the 
retarded T-matrix, we need to evaluate Gq(z) and Fq(z) 



in (21) for z = cj + zO + . Both Green's functions have real 



and imaginary parts: 

G (u>) = G' (u>)+iG%(uj) 

F H = ^H + ^H. 



(22) 



It is convenient to express all energies in units of the 
zero- field magnon bandwidth, Wo = 4 Jo 5, and the field 
in units of the saturation field of the host, h s = SJoS. 
Thus, in the following h = h/h s and Co = uo/AJqS. Three 
other energy scales are needed for the results below: the 
field- dependent gap in the magnon spectrum 



wgap 



2/i, 



an auxiliary scale 



4J S y/i _ 2h? ' 
and the field-dependent magnon bandwidth 



W 
4Jo^ 



„ h < l/y/3 
, h>l/V3. 



(23) 



(24) 



(25) 



For — W < uo < W, and after some algebra, we arrive 
at the following expressions for the imaginary parts 



irWor(uj) 



-q-fe 2 ) 

irWor(Cd) 



Q{oJ 2 -Co 2 g ) (l + oj + h^Ki^) 



+ (l + oj + h 2 l2 )K( 1 ' 2 ) 



(26) 



e(w 2 -^) 7l if( 7 ;) + 72 x( 7 ^) 



The expressions for the real parts of the local Green's 
functions for the same energy range — W < uo < W are 
complimentary to ([26]) 



(20) G' (u) 



e(co 2 -Co 2 g ) (l + LO + h 2 7l )K( 7l ) 



irW r(uj) 

+ e {uj 2 - w 2 ) (i + q + h 2 7l ) —K (— 

y 7i V7i 



sig: 



n( 72 ) (l+6j + h 2 j 2 )K(~f 2 ) 



(27) 



~(l-fe 2 ) r 
7rWor(Cd) 



e(u; 2 -^ 2 ) 7l K( 7l ) 



6 {uo 2 g - uo 2 ) K 



7i 



sign (72)72^(72) 



where O's are step- functions, if s are complete elliptic 
integrals of the first kind, 7^ = yj 1 — 7 2 , and 



7i,2 



h 2 ± r(Cd) 
l-2h 2 ' 



are the roots of the equation uo 2 



0, i.e. 



r{uo) = J(l-2h 2 ){ 



■Co 2 ). 



(28) 



(29) 
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One may notice, that for the field h > 1/y/S some of the 
step-functions [0 (a) 2 — Cj 2 )] are zero in the considered 
energy range uj 2 < W 2 , because W = u gap for this field 
range. As a consequence, some of the terms in (26) and 



le-9 



(27) vanish entirely for that field range. 

For energies outside the magnon bandwidth, uo 2 > W 2 , 
the imaginary parts of Go and Fo are identically zero. 
The derivation of the real parts, while somewhat more 
convoluted, eventually leads to two cases. First, 



G' (w) = 

f oM = 

with 



7rWo^l(^) 

2(1 -h 2 ) 
7rWori(uj) 



Im 



Im 



(i 



i 

7i 



-K 



K 



r 1 (a)) = A /(l-2ft2)(a>2_ fi) 2 i) . 



(30) 



(31) 



This is valid for fields within the range h < l/y/3 for 
any uj 2 

(uj 2 m >W 2 ). Second, 



> W 2 and for 1/^3 < h < l/x/2 for uj 2 > uj 2 



G' (u>) 



i 



7rWor(uj) 



(i 



7i V7i 



(I+W + P72) — K 

72 



(32) 



(1 ~ h 2 ) 
7rWor((d) 



which is valid for fields within the range l/y/3 < h < 1 
for any cj 2 > W 2 and for the fields l/y/3 <h< l/y/2 for 
c^<^ 2 < W 2 . 

Fig. [8] shows the real and imaginary parts of Go and Fq 
for a representative choice of h = Jo S (h = 0.125). Pro- 
nounced van Hove singularities at the top of the magnon 
spectrum (W « 470^) and at the field-induced gap 




FIG. 8. (color online) Real and imaginary parts of Gq(uj) and 
F (uj) for h = J S (h = 0.125). 




FIG. 9. (color online) The relative difference between the 
imaginary part of the analytical (a) and numerical (n) re- 
solvents £ n (z) for the finite N x N cluster with N = 64, 
exchanges J — Jo — 1, spins S = S' = 1, field h = 2, and 
z — ou + ir] with 77 = 0.005. For these parameters, the uniform 
canting angle is = 0.2527 and the impurity canting angle is 
0i 1.122. 



(wgap — h) are clearly visible. These nonanalyticities 
are present in the cj-dependence of the analytical result 
for — Imt 11 (cj) from (19), which is depicted in Fig. 3 of 



the main text, but are much less pronounced. 



Comparison of analytical and numerical results 

The analytical solution of the artificial problem can 
also be used to check the numerical solution beyond the 
discussion in the main text. To this end, rather than 



performing the integration in (21) exactly, the momen- 
tum sum is carried out numerically on a finite cluster of 
the same size as the one used in the numerical proce- 
dure. The resulting two T-matrices can then be com- 
pared. A typical case is shown in Fig. [9j where we 
depict the relative difference of the diagonal elements 
t n (z) of the resolvents obtained from each of the two 
approaches, versus frequency for a finite N x N cluster 
with TV = 64, exchanges J = Jq = 1, spins S = S' = 1, 
and field h = 2. For this plot, the imaginary broaden- 
ing in z — uj + if] has been set to r] = 0.005. For this 
system size and the choice of parameters, such value of 77 
is small enough, so that each individual delta-function in 
the spectrum, corresponding to every single eigenenergy, 
is resolved. Therefore, this figure demonstrates an agree- 
ment between the two approaches not only to within the 
precision of the linear algebra routines that are used for 
the numerical Bogolyubov transformation, but also on a 
level of resolution down to each individual eigenstate. 
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